J. Phys. B: At. Mol. Opt. Phys. 11, (2008) 1-000 
Printed in Taiwan 



1 



Driven Morse Oscillator: 
Model for Multi-photon Dissociation of 
Nitrogen Oxide 

oo 

■ By Julian Juhi-Lian Ting 

Low-Energy Physics Institute, Tsing-Hua University. 
Hsin-Chu, Taiwan 30043, Republic of China [j] 



3 

in 



> 
o 



(Received 5 February 2008) 



Within a one-dimensional semi-classical model with a Morse potential the possibility 
of infrared multi-photon dissociation of vibrationally excited nitrogen oxide was studied. 
The dissociation thresholds of typical driving forces and couplings were found to be 
similar, which indicates that the results were robust to variations of the potential and of 
the definition of dissociation rate. 
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1. Introduction 

^T>' Experimental microwave multi-photon dissociation of nitrogen oxide (NO) is not yet 

I ■ possible, whereas the corresponding optical photodissociation is readily observed via elec- 

tronic transitions. However, it is important to know the possibility and how one can make 
it, firstly because NO is abundant in nature and in living bodies, and secondly because 
the intensity of lasers has reached the order of magnitude needed to do the experiment. 
Although hydrogen fluoride (HF) has a smaller dissociation energy, from an experimental 
view, because of the stability, simplicity and well documented properties of NO, easier 
preparation of intense molecular beams of NO than of HF, and more efficient probing, 
NO has been chosen for investigation instead of HF. Although experiments were formerly 
conducted as NO mostly in the ground state prepared by supersonic expansion, it has be- 
come possible to transfer population efficiently to other states: Yang et al. (1990) report 
experiments in which NO molecules were populated into an initial vibrational state as 
great as n = 25 by stimulated-emission pumping (SEP) whereas Schiemann et al. (1993) 
showed recently how to populate to a state n = 6 by a stimulated Raman process in- 
volving adiabatic passage (STIRAP) method using pulsed lasers. In our study we chose 
parameters for NO whenever molecular properties were needed for our calculations, so 
that future experiments can be easily compared with the results presented here. The 
behaviour of NO is expected to be typical of that of other diatomic molecules. 

Theoretically much previous work concentrated on the corresponding atomic systems. 



Similar to the present work Blumcl and Smilansky (1987) did a calculation on the hydro- 



gen atom and Broeckhove (1992) made a calculation on HF. But no previous calculation 
of the dissociation thresholds of every initial eigenstates has been performed. In this pa- 
per the Morse potential-energy function for the inter-a tomic interact ion was used; only 



the vibrational levels were considered presently because Morse (1929) and Pekeris (194G) 
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parameter symbol value 

range parameter of the potential a 2.7675 x 10 8 cm" 

reduced mass m 7.4643w 

dissociation energy D 6.4968eF 

Morse frequency loq = \j2Do?- jm 1904.2cm -1 

Table 1. Typical parameters for NO. 



have pointed out that the rotational influence on vibration can be taken into account by 
suitably adjusting the Morse parameters. 



The Morse potential was used by Lennard- Jones and Strachan (1935), Strachan (1935) 



to investigate the interaction of atoms and molecules with solid surfaces in the early days. 



Isnor and Richardson (1971), Walker and Preston (1977) studied the driven Morse oscil- 



lator as a model for infrared multi-photon excitation and for dissociation of molecules, 
stimulated by the possibility of laser isotope separation and bond-selective chemistry. 
The Morse potential lacks the complication of infinitely many bound eigenstates and the 
singularity of the Coulomb potential at the origin of the hydrogen atom. 

The methods to calculate multi-photon dissociation are generally divided into non- 
linear classical-mechanical methods, like that of Goggin and Milonni (1988) and semi- 



quantum methods like that of Heather and Mctiu (1987), Heather and Metiu (1988) 



Tanner and Maricq (1988)| , [Tanner and Maricq (1989) 
this study. 



The latter approach is used in 



2. Model 

We consider an isolated non-rotating NO molecule interacting with a plane-polarised 
harmonic laser field. The dimensionless Hamiltonian of a free Morse oscillator for a 
diatomic molecule is 

in which z — a(r — r e ) denotes displacement of inter-atomic distance from equilibrium, 
p = p 2 /v2mD, and p z is the momentum conjugate to z. The parameters of NO are 



given by Huber and Hcrzberg (1979) and summarised in the TABLE. The Hamiltonian 



for the Morse oscillator in the presence of a typical harmonic laser field reads 

JO 

H = H - (2) 

in which fi is the dipole moment operator. In the formula above two dimensionless 
variables are introduced, namely Q — ujl/luq with u>l the laser frequency, and the di- 
mensionless field strength A — gE^/ailD with El the field strength of the laser and q is 
the effective dipole charge of the molecule. In the numerical calculations below various 
types of couplings are considered. 



The discrete eigenfunctions of the free Morse oscillator were given by Morse (1929) 



and denoted by \n > with the corresponding dimensionless eigenvalues equal to 

(n+|) (n + ±f 

En ~ — 2A^' (3) 

in which A = \/2mD/ah = 55.04. The term multi-photon is generally applied to this 
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system, as each photon has energy about 0.24eV near the Morse frequency (u>o), whereas 
the energy difference between the ground state and the first excited state of the free 
Morse oscillator is about 0.23eV; about 28 photons must be absorbed for the transition 
from the ground state to the continuum. 

Our main concern is the transition between each excited state of the free Morse oscil- 
lator and dissociation. With various initial states \n >, driving frequencies and driving 
amplitudes the dissociation periods required might change dramatically. The problems 
of this calculation are whether our model can really represent the molecular system for 
which we intend it for and what is our definition of dissociation. How robust are the 
results? To address these questions we present two typical definitions of dissociation and 
driving forces. 



3. Numerical results 

The following numerical results were obtained using a method described and verified 
previously by Ting et al. (1992), which is a fast Fourier-transformed grid method also 
considered by Feit et al. (1982), Leforestier et al. (1991). 

3.1. Dipolar coupling 



Following Walker and Preston (1977) , Goggin and Milonni (1988) on a driven Morse os- 



cillator, as a consistency check in the first instance, a cosinusoidal wave form was used. 
Hence the Hamiltonian in Eq. (|J) reads 

AO, 

H = H - — zcos(nt). (4) 



According to Goggin and Milonni (1988), Chclkowski and Bandrauk (1990)| for quantum 



chaos the definition of dissociation rate is 

Pdi SS (t) = l-J2\ <r #(*) > I 2 - (5) 

n 

We define Pdissif) — 0.05 for the molecule to satisfy the dissociation condition. If the 
molecule cannot reach such a dissociation rate within 300 cycles we consider it to be 
not dissociable under such conditions of field intensity and frequency. This definition 
is almost equivalent to comparing the dissociation rate after 300 cycles. According to 
Leforestier et al. (1991) the latter should be easier to calculate. However, these two def- 
initions are distinct. Furthermore, the definition of Pdiss is only approximate because, 
when the field is applied, the projection onto the bound states fails to represent the 
true bound-state population unless one assumes that the electro-magnetic field can be 
terminated instantaneously. The values 300 and 0.05 are set according to previous ex- 
perience, although there are some guiding principles: 300 cycles corresponds to a pulse 
of duration about 5 picosecond, which is the time scale for bound dissociation to oc- 
cur. These choices are analogous to those used in the calculation of atomic ionisation by 



Bliimcl and Smilansky (1987). 

The dissociation histories from Eq. (|^) were calculated at fl = 0.9. Fig.l shows disso- 
ciation histories of four types at critical amplitudes, i.e., the largest amplitude at which 
the molecule just fails to be dissociated within 300 vibrations. There are several such 
critical amplitudes for some low- lying initial states ( vide infra ) . This classification is not 
rigorous, but represents at least some typical dissociation histories; it results from inves- 
tigation of more than 500 curves obtained for the calculation for Fig. 3. Ting et al. (1992) 
shows other curves. A brief description of each type follows. For the first type the curve 





Figure 1. Dissociation histories of four types observed at critical amplitudes A — A c and laser 
frequency Q = 0.9; (a) the curves rise smoothly; (b) the curves rise sharply; (c) curves that are 
distinct from those of neighbouring states; (d) straight lines. 



rises smoothly, whereas for the second type the curve rises sharply. For those initial 
states with two transitions the smaller value corresponds to the first type whereas the 
larger value of A c corresponds to the second type. All those states belonging to the 
second type are expected to have a smaller value of A c , although only part of them 
were found. The third type is essentially the same as the first type but their neighbours 
belong to the second type. The fourth type belongs to the initial states exceeding 44; 
the dissociation history is almost a straight line. Except for the latter type, for which 
the figure cannot display clearly, the reason for a greater dissociation period is not only 
that the slope of the curve was altered, but that there were several flat regions in which 
there was little molecular dissociation. As a further investigation for the third type, Fig. 2 
compared at the same amplitude with initial states near no = 28, indicating that some 
initial states such as no = 28 dissociate more readily than their neighbours. These states 



correspond to the ionisation windows proposed by Blumel and Smilansky (1987). There 



are indications that some states dissociate more easily than others. The dissociation 
process might be taken to have two steps; i.e., the population first transferred to those 
easy states then to other states. Generally the populations deployed themselves soon af- 
ter the process began and stayed well localised, as Walker and Preston (1977) observed. 
3hankar (1980), Pfcifcr (1993) discussed that an initial set-up time is required for any 
periodic system is a general property, because the system has to sense the periodic per- 
turbation. Various quantities such as the inter-atomic distance of initial states and the 
average variation of energy with time were also computed. A preliminary result on the 
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Amplitude 



Fig. 3 




Initial State 



Figure 3. Critical amplitude {A c ) for various initial states at Q — 0.9; the Hamiltonian is Eq. 
(0) and the dissociation rate is Eq. (|B|); for low-lying states which have two values of A c only 
the lower one is connected. 



inter-atomic distance shows that, immediately after the field is switched on, states with 
no > 25 move in a direction distinct from that for no < 25 away from the inter-atomic 
equilibrium position. 

In Fig. 3 the critical fields required for the molecule to become dissociable is plotted. 
For multi-critical states the lowest ones are connected and the others marked. Therefore 
the plot provides only a sufficient criterion for dissociation, not a necessary condition. 
The reason is that as we use numerical methods we can know only the points that we 
computed. The amplitude resolution of our calculation is AA = 0.0005. If there were 
a narrow transition less than 0.0005 one might not be able to find it, because such a 
point would imply experimentally a laser of precise power. We find at most two critical 
amplitudes in the present case while in another study using two-color laser fields we find 
up to three. The multi-critical amplitude is also a result of our selection of the critical 
dissociation period to be 300. However, this point implies that decreasing field strength 
does not necessarily increasing the dissociation period. Lu et al. (1991) noted that the 
finite size of H probably places a characteristic scale in parameter space less than which 
the dissociation boundary ceases to show fractal structure. Another feature of the figure 
is that the largest value of A c is several orders of magnitude larger than the smallest one. 
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Figure 4 
d (b) Q 



. Critical amplitude (A c ) for various initial states at (a) £1 — 0.89 (solid thick line) 
= 0.91 (thin dashed line); the Hamiltonian is Eq. (||) and the dissociation rate is Eq. 
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Figure 5. Critical amplitude (A c ) vs. frequency of initial states (a) no = 1, (b) no = 17 and 

(c) n = 24. 



A calculation shows that for no > 6 the photon energy becomes larger than the energy 
difference between adjacent states. 

To discover whether = 0.9 is the most efficient frequency for all initial states we 
tested ft = 0.89 and £1 — 0.91. The critical amplitudes are plotted in Fig. 4. As we had 
insufficient computer time to test all states and all frequencies, we selected n Q = 1, 17 
and 24 to find whether values of A c alter for frequencies from 0.8 to 1.1 at interval 0.01. 
The plots appear in Fig. 5. 

3.2. Exponential coupling 

The perturbation in Eq. ([|), although simple and easily done by the Floquet method, has 
several problems in relation to reality. Firstly, the molecular dipole and laser coupling in 
equation Eq. (|J) assumes a infinite range of constant coupling, which is satisfactory for 
an atomic system but not for a real molecule. Secondly, as Sakurai (1967)| has noticed, 



the field should be zero at t = and gradually increases from zero to a finite value; this 
problem is called adiabatic switching. Therefore, another exponential form of coupling 

H = H ^(z + a)e-^/ b sin(ffi) sin 2 (^), (6) 
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Figure 6. Critical amplitude (A c ) of various initial states and modified potential Eq. (g) at 
Q = 0.9; for the solid thick line T — 300 whereas for the thin dashed line T = 600. 



was used to test whether there is any effect of the approximation. The problem of the 
definition of dissociation rate was point ed out above (below Eq.(||)); here we alter it to 
comply with Heather and Mctiu (1987)| 

Pdissit) = 1- <if)(t)\i/)(t) > . (7) 
This definition is appropriate if coherent excitation becomes important. The variables 



and b were chosen to have t h e values 2 and 1 respectively in the calculation. Heather and Mctiu (1987)| . 



Heather and Metiu (1 



Tanner and Maricq (1 



Tanner and Maricq (1989) have 



also used coupling functions of similar form. We take T, the pulse duration, to be 300 
optical cycles. This choice has the advantage that, when A = A e , Td is almost 300 cycles 
and the field strength is almost zero at that time. Hence there is little effect of the defini- 
tion of dissociation rate. The resulting critical front appears in Fig. 6 as a solid thick line. 
There is no essential distinction from Fig. 3 and Fig. 4 except a scaling parameter that 
can be adjusted by choosing a and b. For lower initial states there is more variation than 
for higher initial states. To assess whether T is large enough, a longer period T = 600 is 
tested. The result is plotted in Fig. 6 with a thin dashed line; the two lines show similar 
behaviour. 



4. Summary 

In conclusion, thresholds of vibrational dissociation of NO are considered in this paper. 
The main result of this paper is depicted in Fig. 6. Two perturbations and definitions of 
dissociation rate lead to similar curves with the transition from high A c to low A c shifted 
slightly to higher initial states. Laser fields and couplings of other kinds were tested for 
which similar curves were obtained. Therefore, as our main concern is the critical ampli- 
tude which is an indirect result, we conclude that our results are robust to variation of 
the potential function. Variations of other types such as the laser duration and definition 
of dissociation threshold are less significant. A rough estimate of the decrease of A c is 
given. The assumption behind this model is that no other electronically nonadiabatic 
effects will occur first. One must take care in comparing the results presented in this 
work with those from experiment. For instance, because of the finite cross section of the 
beam, not all molecules experience the same field strength. The ignorance of rotational 
influence should also be remembered. The additional selection rule in the rotating case 
may make the progression to dissociation harder. 
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